The build-up, configuration, and dynamical sensitivity of the Eurasian ice-sheet complex to Late Weichselian climatic and oceanic forcing

The Eurasian ice-sheet complex (EISC) was the third largest ice mass during the Last Glacial Maximum (LGM), after the Antarctic and North American ice sheets. Despite its global significance, a comprehensive account of its evolution from independent nucleation centres to its maximum extent is conspicuously lacking. Here, a first-order, thermomechanical model, robustly constrained by empirical evidence, is used to investigate the dynamics of the EISC throughout its build-up to its maximum configuration. The ice flow model is coupled to a reference climate and applied at 10 km spatial resolution across a domain that includes the three main spreading centres of the Celtic, Fennoscandian and Barents Sea ice sheets. The model is forced with the NGRIP palaeo-isotope curve from 37 ka BP onwards and model skill is assessed against collated flowsets, marginal moraines, exposure ages and relative sea-level history. The evolution of the EISC to its LGM configuration was complex and asynchronous; the western, maritime margins of the Fennoscandian and Celtic ice sheets responded rapidly and advanced across their continental shelves by 29 ka BP, yet the maximum aerial extent (5.48 × 106 km2) and volume (7.18 × 106 km3) of the ice complex was attained some 6 ka later at c. 22.7 ka BP. This maximum stand was short-lived as the North Sea and Atlantic margins were already in retreat whilst eastern margins were still advancing up until c. 20 ka BP. High rates of basal erosion are modelled beneath ice streams and outlet glaciers draining the Celtic and Fennoscandian ice sheets with extensive preservation elsewhere due to frozen subglacial conditions, including much of the Barents and Kara seas. Here, and elsewhere across the Norwegian shelf and North Sea, high pressure subglacial conditions would have promoted localised gas hydrate formation

Nunataks as barriers to ice flow : implications for palaeo ice sheet reconstructions

Funding: This research has been supported by the Vetenskapsrådet (grant no. 2016-04422), the Deutsche Forschungsgemeinschaft (grant no. 1158-365737614), the National Science Foundation (grant no. OPP-1542930), and the Norsk Polarinstitutt (grant no. 2015/38/7/NK/ihs).Numerical models predict that discharge from the polar ice sheets will become the largest contributor to sea-level rise over the coming centuries. However, the predicted amount of ice discharge and associated thinning depends on how well ice sheet models reproduce glaciological processes, such as ice flow in regions of large topographic relief, where ice flows around bedrock summits (i.e. nunataks) and through outlet glaciers. The ability of ice sheet models to capture long-term ice loss is best tested by comparing model simulations against geological data. A benchmark for such models is ice surface elevation change, which has been constrained empirically at nunataks and along margins of outlet glaciers using cosmogenic exposure dating. However, the usefulness of this approach in quantifying ice sheet thinning relies on how well such records represent changes in regional ice surface elevation. Here we examine how ice surface elevations respond to the presence of strong topographic relief that acts as an obstacle by modelling ice flow around and between idealised nunataks during periods of imposed ice sheet thinning. We find that, for realistic Antarctic conditions, a single nunatak can exert an impact on ice thickness over 20 km away from its summit, with its most prominent effect being a local increase (decrease) of the ice surface elevation of hundreds of metres upstream (downstream) of the obstacle. A direct consequence of this differential surface response for cosmogenic exposure dating is a delay in the time of bedrock exposure upstream relative to downstream of a nunatak summit. A nunatak elongated transversely to ice flow is able to increase ice retention and therefore impose steeper ice surface gradients, while efficient ice drainage through outlet glaciers produces gentler gradients. Such differences, however, are not typically captured by continent-wide ice sheet models due to their coarse grid resolutions. Their inability to capture site-specific surface elevation changes appears to be a key reason for the observed mismatches between the timing of ice-free conditions from cosmogenic exposure dating and model simulations. We conclude that a model grid refinement over complex topography and information about sample position relative to ice flow near the nunatak are necessary to improve data–model comparisons of ice surface elevation and therefore the ability of models to simulate ice discharge in regions of large topographic relief.Publisher PDFPeer reviewe

Deglaciation of Fennoscandia

To provide a new reconstruction of the deglaciation of the Fennoscandian Ice Sheet, in the form of calendar-year time-slices, which are particularly useful for ice sheet modelling, we have compiled and synthesized published geomorphological data for eskers, ice-marginal formations, lineations, marginal meltwater channels, striae, ice-dammed lakes, and geochronological data from radiocarbon, varve, optically-stimulated luminescence, and cosmogenic nuclide dating. This 25 is summarized as a deglaciation map of the Fennoscandian Ice Sheet with isochrons marking every 1000 years between 22 and 13 cal kyr BP and every hundred years between 11.6 and final ice decay after 9.7 cal kyr BP. Deglaciation patterns vary across the Fennoscandian Ice Sheet domain, reflecting differences in climatic and geomorphic settings as well as ice sheet basal thermal conditions and terrestrial versus marine margins. For example, the ice sheet margin in the high-precipitation coastal setting of the western sector responded sensitively to climatic variations leaving a detailed record of prominent moraines and ice-marginal deposits in many fjords and coastal valleys. Retreat rates across the southern sector differed between slow retreat of the terrestrial margin in western and southern Sweden and rapid retreat of the calving ice margin in the Baltic Basin. Our reconstruction is consistent with much of the published research. However, the synthesis of a large amount of existing and new data support refined reconstructions in some areas. For example, we locate the LGM extent of the ice sheet in northwestern Russia further east than previously suggested and conclude that it occurred at a later time than the rest of the ice sheet, at around 17-15 cal kyr BP, and propose a slightly different chronology of moraine formation over southern Sweden based on improved correlations of moraine segments using new LiDAR data and tying the timing of moraine formation to Greenland ice core cold stages. Retreat rates vary by as much as an order of magnitude in different sectors of the ice sheet, with the lowest rates on the high-elevation and maritime Norwegian margin. Retreat rates compared to the climatic information provided by the Greenland ice core record show a general correspondence between retreat rate and climatic forcing, although a close match between retreat rate and climate is unlikely because of other controls, such as topography and marine versus terrestrial margins. Overall, the time slice reconstructions of Fennoscandian Ice Sheet deglaciation from 22 to 9.7 cal kyr BP provide an important dataset for understanding the contexts that underpin spatial and temporal patterns in retreat of the Fennoscandian Ice Sheet, and are an important resource for testing and refining ice sheet models

Changes in vertical ice extent along the East Antarctic Ice Sheet margin in western Dronning Maud Land – initial field and modelling results of the MAGIC-DML collaboration

Constraining numerical ice sheet models by comparison with observational data is crucial to address the interactions between cryosphere and climate at a wide range of scales. Such models are tested and refined by comparing model predictions of past ice geometries with field-based reconstructions from geological, geomorphological, and ice core data. However, for the East Antarctic Ice sheet, there is a critical gap in the empirical data necessary to reconstruct changes in ice sheet geometry in the Dronning Maud Land (DML) region. In addition, there is poor control on the regional climate history of the ice sheet margin, because ice-core locations, where detailed reconstructions of climate history exist, are located on high inland domes. This leaves numerical models ofregional glaciation history in this near-coastal area largely unconstrained. MAGIC-DML is an ongoing Swedish-US-Norwegian-German-UK collaboration with a focus on improvingice sheet models of the western DML margin by combining advances in modeling with filling critical data gaps regarding the timing and pattern of ice-surface changes. A combination of geomorphological mapping using remote sensing data, field observations, cosmogenic nuclide surface exposure dating, and numerical ice sheetmodeling are being used in an iterative manner to produce a comprehensive reconstruction of the glacial historyof western DML. Here, we present an overview of the project, field evidence for formerly higher ice surfaces and in-situ cosmogenic nuclide measurements from the 2016/17 expedition. Preliminary field evidence indicate that interior sectors of DML have experienced a general decrease in ice sheet thickness since the late Miocene, with potential episodes of increasing thickness in the late Pleistocene (700-300 ka, 250-75 ka). To aid in interpreting these field data, new high-resolution ice sheet model reconstructions, constraining ice sheet configurations during key episodes, are presented

Mid-Pleistocene ice sheet fluctuations from cosmogenic nuclide field constraints in western Dronning Maud Land, Antarctica

The East Antarctic Ice Sheet (EAIS) is generally assumed to have been relatively insensitive to Quaternary climat echange. However, recent studies have shown potential instabilities in coastal, marine sectors of the EAIS. In addition, long-term climate reconstructions and modelling experiments indicate the potential for significant changes in ice volume and ice sheet configuration since the Pliocene. Hence, more empirical evidence for ice surface and ice volume changes is required to discriminate between contrasting inferences. MAGIC-DML is an ongoing Swedish-US-Norwegian-German-UK collaboration focused on improving ice sheetm odels by filling critical data gaps that exist in our knowledge of the timing and pattern of ice surface changes along the western Dronning Maud Land (DML) margin and combining this with advances in numerical techniques. As part of the project, field studies in the 2016/17 and 2017/18 austral summers targeted selected sites spanning accessible altitudes in the Heimefrontfjella, Vestfjella, Ahlmannryggen, Borgmassivet, and Kirwanveggen nunatakranges for in situcosmogenic nuclide sampling. Comparing concentrations of nuclides with widely differing half-lives in bedrock and erratics from a range of altitudes above modern ice surfaces can provide information on ice sheet fluctuations and complex burial and exposure histories, and thus, past configurations of non-erosive ice. Quartz-bearing rock types were sampled and analyzed for 10Be (t1/21.4 My),14C (t1/25.7 ky),26Al (t1/2705ky), and 21Ne (stable), and mafic lithologies for36Cl (t1/2301 ky). Results thus far for 3210Be and 26Al isotope pairs complemented with seven21Ne measurements have yielded some consistent patterns of paleoglaciation for the western DML margin. Eight out of fourteen bedrock samples from high-elevation (1700-2238 m a.s.l.) ridges and summits return some of the oldest exposure ages in Antarctica and have consistent 10Be,26Al, and 21Ne minimum apparent exposure ages of 1.8-4.1 Ma. Initial results therefore indicate that parts of the ice sheet marginal to the Antarctic plateau, along the Heimefrontfjella range, generally have experienced a decrease in ice thickness since the late Miocene. Another six bedrock samples (1556-1732 ma.s.l.) fall in the 300-700 ka range, and they all show significant burial. At face value, perhaps this indicates aregional ice sheet surface above 1700 m a.s.l. for much of the Plio-early Pleistocene. All other samples analyzedto date are erratics from lower elevation and more coastal sites (10 from nunataks at 553-1400 m a.s.l., and 6 froma surface moraine at 1385 m a.s.l.), exhibiting ages between 59 and 275 ka, save for two (4 and 6 ka). Whereas almost all of the nunatak erratics (including the young ones) show significant burial durations, five of the six surface moraine samples do not. These 2016/17 field samples are not yet leading to conclusive age constraints but already start to paint a picture of the western DML margin being relatively stable although there was possibly one or more episodes of relatively limited ice thickening during the last 700 ka

Glacial geomorphology of the Bayan Har sector of the NE Tibetan Plateau

We here present a detailed glacial geomorphological map covering 136,500 km2 of the Bayan Har sector of the northeastern Tibetan Plateau - an area previously suggested to have nourished the most extensive Quaternary glaciers of the Tibetan Plateau. The map, presented at a scale of 1:650,000, is based on remote sensing of a 90 m SRTM digital elevation model and 15/30 m Landsat ETM+ satellite imagery. Seven landform types have been mapped; glacial valleys, glacial troughs, glacial lineations,marginal moraines, marginal moraine remnants, meltwater channels and hummocky terrain. A large number of glacial landforms exist, concentrated around mountain blocks protruding above the surrounding plateau area, testifying to former glacial activity. In contrast, large plateau areas of lower altitude lack glacial landforms. The mapped glacial geomorphology indicates multiple former glacial advances primarily by valley and piedmont glaciers, but lends no support to the hypothesis of ice sheet scale glaciation in the area. The presented glacial geomorphological map demonstrates the usefulness of remote sensing techniques for mapping the glacial geomorphology of the Tibetan Plateau, and it will be used for reconstructing the paleoglaciology of the Bayan Har sector of the northeastern Tibetan Plateau.Part of urn:nbn:se:su:diva-750

Inceptions: mechanisms, patterns and timing of ice sheet inception

No abstract available

Glacial geomorphology of the central Tibetan Plateau

The glacial geomorphology of the central Tibetan Plateau was mapped over 285,000 km2. Here we present a map covering 135,000 km2 at a scale of 1:660,000. The glacial geomorphology was mapped using 15 and 30 m resolution Landsat 7 ETM+ satellite imagery, a 90 m resolution SRTM digital elevation model, and satellite and aerial images displayed in Google Earth. Four landform categories were discernible and mapped; glacial valleys, marginal moraines, glacial lineations, and hummocky terrain. The mapped landforms indicate multiple glacial advances of valley and piedmont glaciers. The mapped landform record lends no support to individual ice centres having coalesced to form a plateau-wide ice sheet.